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Abstract 

We compared the brain areas that showeci significant flow changes induced by 
selective stimulation of muscle and cutaneous afferents using fMRI BOLD 
imaging. Afferents arising from the right hand were studied in eight volun- 
teers with electrical stimulation of the digital nerve of the index finger and 
over the motor point of the FDI muscle. Both methods evoked areas of signif- 
icant activation corticaUy, subcortically, and in the cerebellum. Selective mus- 
cle afferent stimulation caused significant activation in motor-related areas. It 
also caused significantly greater activation within the contralateral precentral 
gyrus, insula, and within the ipsilateral cerebellum as well as greater areas of 
reduced blood flow when compared to the cutaneous stimuli. We demon- 
strated separate precentral and postcentral foci of excitation with muscle affer- 
ent stimulation. We conclude, contrary to the findings with evoked potentials, 
that muscle afferents evoke more widespread cortical, subcortical, and cerebel- 
lar activation than do cutaneous afferents. This emphasizes the importance, 
for studies of movement, of matching the kinematic aspects in order to avoid 
the results being confounded by alterations in muscle afferent activation. The 
findings are consistent with clinical observations of the movement 
consequences of sensory loss and may also be the basis for the contribution of 
disturbed sensorimotor processing to disorders of movement. 
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Introduction 

Both muscle and skin are richly innervated with 
specialized receptors which provide sensory information to 
the central nervous system. While many nerves are mixed, 
there exist both primary sensory and primary motor 
(including afferents) nerves. One-third to one-half of the 
myelinated fibers in mammalian muscle nerves are afferent 



(Sherrington 1894). Projections to the somatosensory 
cortex have been previously studied using direct electrical 
stimulation of peripheral nerves. Electrical stimulation of 
the median nerve depolarizes motor axons and large diam- 
eter, fast conducting sensory fibers of types Aa and Afi 
(Groups I and II). These derive from cutaneous receptors, 
muscle spindles (1° and 2°), Golgi tendon organs, and 
joint capsule receptors (Marani and Lakke 2012). Most 
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Afferent Activation of Sensorimotor Cortex and Cerebellum 

imaging studies have concentrated on the cortical effects, 
particularly for the sensorimotor cortex. Electrical stimula- 
tion of the median nerve, or the application of air puff 
stimuli to the fingers, sequentially activates contralateral SI 
followed by posterior parietal cortex and bilateral SII, as 
measured with magnetoencephalography (Forss et al. 
1994a,b). Functional MRI (1.5 T) has been used with elec- 
trical stimulation of the right median nerve, and has 
shown cortical activity in the SI hand area, with SII being 
activated at lower current values (Backes et al. 2000). 
Krause et al. (2001) in an fMRI (1.5 T) study showed that 
electrical stimuli applied to the fingers thus activating cuta- 
neous nerve afferents only, could selectively activate areas 
3b and the crown of the postcentral gyrus (area 1-2). It is 
also well known that muscle and skin afferent projections 
are not limited to the primary sensorimotor cortex. For 
example, passive joint movements, potentially activating 
both muscle and cutaneous afferents, have been shown to 
excite neurons within the motor cortex, thalamus, and 
basal ganglia (e.g.. Lemon and Porter 1976; Horne and 
Porter 1980; Colebatch et al. 1990; Weiller et al. 1996; 
Lozano and Hutchison 2002), and SII (the secondary sen- 
sory area) is an additional projection area for both muscle 
and cutaneous afferents (Landgren et al. 1967). Joint affer- 
ents are thought to make only a minor contribution as 
they usually discharge only at the extremes of joint move- 
ment (McCloskey 1978). 

The functional roles of cutaneous and muscle receptors 
are different. Low threshold cutaneous receptors for 
example, can provide information about the light touch 
of a caress, the texture of a surface and the presence and 
location of an insect crawling on the skin. Muscle 
receptors have a more direct role in motor control - for 
example in the production of a constant force contraction 
(Rothwell et al. 1982; Proske and Gandevia 2012). Both 
afferent species have a role in proprioception, particularly 
for the hand (Collins et al. 2000; for review Proske and 
Gandevia 2012). This important sense allows an individ- 
ual to recognize the position and movement of the limbs 
and is required for the production of accurate and finely 
graded contractions. Our study was designed to compare 
the patterns of cortical activation evoked by a cutaneous 
versus a muscle afferent volley arising from the same 
region of the hand to examine whether these functional 
differences are reflected in the patterns of evoked suprasp- 
inal excitation. 

Methods 

Subjects 

Eight male adult volunteers were recruited from the 
personnel employed at, or associates of, the then Prince of 
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Wales Medical Research Institute (now Neuroscience 
Research Australia). The mean age (±SD) of the subjects 
was 40.4 ± 8.8 years (27-52 years) and the mean weight 
79.3 ± 12.1 kg (55-95 kg). One subject was left-handed. 
AU subjects had no past medical history of neurological or 
psychiatric disease. All subjects signed a written statement 
of informed consent approved by the University of New 
South Wales Human Ethics Research Committee (HREC 
03255). 

Stimulus paradigm 

Electrical stimulation was applied via carbon electrodes 
(~3 cm^ Red Dot™ 3M, St Paul, MN) at two sites on the 
right hand. The anode for both cutaneous and muscle 
stimuli was applied to the lateral skin (thumb side) of the 
right index finger approximately 1-2 cm distal to the 
metacarpophalangeal joint (Fig. 1). The cutaneous stimu- 
lus cathode was applied on the medial side of the index 
finger 1-2 cm distal to the metacarpophalangeal joint. 
The muscle stimulus cathode was applied over the motor 
point of the first dorsal interosseous (FDI) muscle after 
delineation with a probe electrode. The skin over the 
motor point of FDI was anesthetized using a topical local 
anesthetic (EMLA® AstraZeneca, Sydney, Australia: 
Lignocaine 25 mg/g and Prilocaine 25 mg/g) applied 
under nonporous occlusion (Tegaderm™ 3M, St Paul, 
MN) for 30-45 min prior to electrode placement to avoid 
direct cutaneous afferent stimulation. AU electrodes were 
secured with adhesive tape (Micropore™ 3M, St Paul, 
MN). The electrical stimulus consisted of rectangular 
pulses at 4 Hz and 0.1 ms duration. Electrical stimulation 
of the finger at frequencies of 4-16 Hz produces good 
activation of the somatosensory cortex (Takanashi et al. 
2003). Stimuli were generated by an isolated current stim- 
ulator (DS7 Digittmer Ltd, UK) located outside the room 
that was triggered by a purpose-built timing circuit at -3/ 
sec. Electrical stimuli were delivered via 5-m twisted cop- 
per cables connected to the carbon electrodes by clips that 
were also secured with adhesive tape (Fig. 1). 

Individual subject thresholds for stimulus current were 
defined as the mean perceptual cutaneous threshold or 
the mean muscle twitch (motor) threshold. The mean 
cutaneous threshold was 7.3 ± 1.5 mA (range 5.5- 
10.0 mA) and the mean muscle stimulus threshold was 
9.3 ± 2.5 mA (range 5.0-12.5 mA). Two different 
stimulus amplitudes were delivered in both the cutaneous 
and muscle conditions: twice (2T) and four (4T) times 
the stimulus threshold. Most subjects perceived the 4T 
muscle stimuli as somewhat painful. 

The experimental design for each subject consisted of 
eight separate fMRI experimental runs as well as both T2 
and 3D volume image scans. The stimulus runs consisted 



2014 I Vol. 2 I Iss. 4 I e00270 
Page 2 



© 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of 
the American Physiological Society and The Physiological Society. 



D. L. Wardman ef al. 



Afferent Activation of Sensorimotor Cortex and Cerebellum 




I 
I 



0 50 100 150 

Tim© (sec) 

Figure 1. Experimental setup and stimulus block design. (Upper 
Inalf) The cutaneous stimulus cathode was applied on the medial 
side of the index finger 1-2 cm distal to the metacarpophalangeal 
joint (black and red electrodes). The muscle stimulus cathode was 
applied over the motor point of the first dorsal interosseus muscle 
after delineation with a probe electrode (circle with black center). 
The stimulus anode in both cases was applied to the lateral skin 
(thumb side) of the right index finger approximately 1-2 cm distal 
to the metacarpophalangeal joint. (Lower half) The trial block 
design consisted of 8 fMRI collection runs. The order of the 
cutaneous or muscle stimulus blocks as well as the order of 
stimulus intensity (2T or 4T) were randomized. The different 
stimulus paradigms started with either a rest period followed by 
stimulation (A) or an initial stimulus (B) period (see text). 

of four fMRI collection periods for both cutaneous and 
muscle stimuli, at either the 2T or 4T stimulus intensity 
followed by a rest period, and with each stimulus 
intensity presented twice. For each specific stimulus, a 
boxcar paradigm was used, consisting of 60 dynamic 
scans (InterScan interval 3 sec), divided into six periods, 
three stimulus periods and three rest periods, with each 
period having a duration of 30 sec (TR 3000 per dynamic 
scan). Two different paradigms, A and B, were used. 



which corresponded to the fMRI collection run 
commencing with either (A) a rest period (with stimulus 
onsets at scans 11, 31 and 51) or (B) a stimulus period 
(stimulus onsets at scans 01, 21 and 41). One of four 
different protocols was randomly used: AAAA; BBBB; 
ABAB; BABA. 

Functional MRI experimental data 
acquisition 

All experimental data were acquired using a 3T MR 
Achieva scanner (Philips) which was equipped with a 
transmit-receive head coil. Image artifacts were avoided 
by stretching the copper stimulation cables as much as 
possible and by using (nonmagnetic) carbon electrodes. 
Subjects were supine and head fixation was achieved by 
two foam cushions and by straps across the head. Subjects 
were instructed to keep their arms relaxed by their sides. 
The scanner room was dimly lit and subjects were 
instructed to relax and close their eyes during the 
experiment but were not permitted to sleep. Subjects 
wore hearing protection and were connected via intercom 
with the control room. The entire study lasted at most 
90 min per subject, including setup time. 

A T2-weighted multishot (TSE) scan preceded the 
functional scans with parameters Tr 4101.0 ms, T^ 
120.0 ms, flip angle 90.0°, field of view 230.0 mm, matrix 
dimensions 256 x 256, 30 contiguous transverse slices and 
a reconstructed voxel size of 0.45 x 0.45 x 4.0 mm. An 
approximately 12-cm thick stack of slices was defined, 
encompassing transverse slices that covered both the cere- 
brum and cerebellum. The functional scan consisted of a 
single shot multiple slice T2* -sensitive echo planar imaging 
(EPI) sequence, with parameters Tr 3000 ms, Te 33 ms, 
flip angle 90.0°, field of view 230.0 mm, matrix dimension 
96 x 128, 30 contiguous slices and a reconstructed voxel 
size of 1.8 X 1.8 x 4.0 mm. The fMRI signal was based on 
the blood-oxygen-level-dependent (BOLD) effect. The 
fMRI sequence incorporated dummy pulses prior to image 
acquisition to ensure adequate steady-state imaging. A 3D 
Tj-weighted multishot (TEE) echo scan succeeded the 
functional scans with parameters Tr 6.7 ms, Tg 3.0 ms, flip 
angle 8.0°, field of view 230.0 mm, matrix dimensions 
256 X 256, 200 contiguous coronal slices, and a recon- 
structed voxel size of 0.9 x 0.9 x 0.9 mm. 

Data analysis 

The acquired image data were converted to Analyze 7.5 
format using a conversion program, MRIconvert program 
(Version 1.02W: J.Smith, Lewis Centre for Neurolmaging, 
University of Oregon). All image data were analyzed using 
the Statistical Parametric Mapping (SPM2: Wellcome 
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Department of Imaging Neuroscience, Institute of Neu- 
rology, University College London) software package 
implemented in MATLAB 6.5.1 (The Mathworks Inc., 
Natick, MA). 

Post processing steps included image realignment, nor- 
malization, and smoothing. The scans were first corrected 
for small motion artifacts by trilinear realignment (quality 
0.75) to the first volume scan of each time series. Subse- 
quently, the images were normalized to the fourth image by 
transformation into standard space, using the EPI template 
image of the Montreal Neurological Institute (MNI). The 
voxel sizes of the normalized images were 2x2x4 mm 
isotropically with trilinear interpolation and template 
bounding box ( - 90:90, - 1 26:90, - 72: 108) . To enhance the 
spatial signal-to-noise ratio and to facilitate intersubject 
averaging, the normalized images were smoothed with a 
Gaussian filter with a fuU width at half maximum of 8 mm. 
Four of the eight runs from one subject were excluded due 
to data corruption during acquisition. 

Data were analyzed for each condition at both the indi- 
vidual and group level as a General Linear Model for fMRI 
time series. The fMRI design matrix for each run was speci- 
fied in scans and vector of onsets, using the hemodynamic 
response function. Runs were not replications due to the 
different onsets for each subject. Smoothed data from each 
run were then scaled before the SPM was estimated. Ini- 
tially a fixed effects analysis was used but this analysis 
showed wide interindividual variation for areas of signifi- 
cant sensorimotor activation. A second-level random 
effects analysis was, therefore, performed. This technique 
allows inferences to be made about the population from 
which the subjects are drawn, rather than just for the group 
of subjects studied (Penny and Holmes 2007). Because we 
were primarily interested in the effects of different types of 
afferents, the main analysis consisted of considering the 
effects of combined muscle (both 2T and 4T) and skin 
(both 2T and 4T) stimulation for all subjects. For the main 
analysis a P < 0.001 (t = 3.73-3.79) criterion was used, 
with a minimum cluster size of 5. The active conditions 
were compared using a random effects analysis using t-tests 
with a less stringent criterion (P < 0.005) but masked by 
the previously demonstrated significant changes for the 
condition of interest. The resulting sets of images repre- 
sented statistical parametric maps of the t statistic SPM{f}. 

The sterotaxic coordinates of Talairach and Tournoux 
(1988) were used to define the Brodmann areas of the 
observed activation foci after transformation from MNI 
coordinates using the Matlab file (mni2tal.m, M. Brett 
1999: Cambridge Imagers Unit, Cambridge, UK). The cor- 
tical areas corresponding to the observed activation foci 
were derived from the Talairach daemon client program 
(Version 2.0; Research Imaging Centre, UTHSCSA). The 
exact or nearest gray matter structure was reported as the 



area corresponding to each coordinate. One point initially 
allocated to area 4 was reassigned to area 3 and one point 
the reverse after review using the Talairach and Tournoux 
(1988) atlas and plotting on a representative high-resolu- 
tion MRI image. All images were displayed according to 
radiological convention, that is the right hemisphere is 
shown on the left side of the image. All tabulated data are 
shown in MNI coordinates. Cerebellar locations were based 
upon the atlas of Schmahmann et al. (1999) and for the 
cerebellar nuclei the work of Dimitrova et al. (2002, 2006). 

Results 

Muscle stimulation (2T and 4T) 

The random effects analysis produced multiple foci of 
significant activation ipsilateraUy in the cerebellum and 
predominantly contralateraUy for the cerebrum and sub- 
cortical areas (Fig. 2, Table 1). Highly significant activa- 
tions were present in the contralateral insula and nearby 
temporal lobe. Activation sites were present within the 
contralateral sensorimotor cortex, in two distinct areas 
corresponding to heights (Z) of 60 and 44-48 mm above 
the reference plane. Activations were also present within 
the contralateral superior temporal gyrus, thalamus, and 
basal ganglia. The significant activations within the 
cerebellum consisted of a near-midline group and a sec- 
ond group lying more laterally deep within the ipsilateral 
cerebellar hemisphere (Table 1). 

The muscle afferent stimulation condition also evoked 
significant areas of reduced flow in a large number of loca- 
tions bilaterally, mainly lying within the posterior parietal 
and occipital lobes and anterior temporal lobes (Table 2). 

Cutaneous stimulation (2T and 4T) 

Cutaneous stimulation activated a wide range of areas, pre- 
dominantly in the contralateral hemisphere (Fig. 3, 
Table 3). Multiple areas of activation within the frontal 
lobes were present, including one area of activation present 
in the contralateral sensorimotor cortex at height 44 mm. 
Although initially allocated by the Talairach daemon client 
program to area 4, review suggested this point was post 
central, on the crown of the gyrus within area 3. Contralat- 
eral parietal lobe activation corresponding to the second 
somatosensory cortex (SII: area 40) was present and there 
was activation within the insula (area 13). Bilateral areas of 
activation were present in the temporal lobes. The cerebel- 
lum showed an area of activation ipsilateraUy with other 
weaker areas of activation, below our statistical criterion. 

Areas of significantly reduced flow were also present 
(Table 4), but substantially fewer than for muscle stimu- 
lation. 
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Z = 48 



Z = 60 




Figure 2. Activations in response to muscle stimulation. Group brain activation in response to stimulation of the motor point of the right first 
dorsal interosseus muscle at 2 and 4T, using random effects analysis. Statistical images were coregistered with the SPM-MNI single subject T1 
image. The axial/horizontal sections in the upper row show activation areas at the sensorimotor cortex. Sections in the middle row show 
activation areas in the second sensory area, thalamus, and insula. Sections in the lower row show activation areas in the right and midline 
cerebellum, including the nucleus (dentate). Colored bar labels indicate levels of activation in t value. For illustration purposes a cutoff of 
f = 3.5 has been used. 



Inter-stimulus comparisons 

Muscle versus cutaneous 

The higher site of sensorimotor cortex activation, 
activation in the contralateral insula, and activation in the 
ipsilateral cerebellum were significantly greater for muscle 
nerve stimulation than for cutaneous nerve stimulation 



(shown in bold in Table 1). Cutaneous stimulation was 
associated with several areas of significantly increased 
activation compared to muscle afferent stimulation, 
including one site within the contralateral superior 
parietal lobule. The other sites lay within the ipsilateral 
hemisphere and included the ipsilateral precentral gyrus 
and multiple sites within the middle temporal gyrus 
(Table 3). 
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Table 2. Muscle stimulation - random effects. Brain areas witli reduced flow in response to stimulation of the rigint first dorsal interosseous 
muscle nerve (t-test, P < 0.001 Uncorrected, Voxels >B). 
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Combined data for both muscle nen/e stimulation intensities - reductions. Threshold for statistical significance was P < 0.001 (uncorrected) 
with a cluster size of 5 or more voxels, using second level random effects analysis. Number of voxels indicates the overall spatial extent of 
significant clusters within the activated volume. Cluster sizes with asterisk are foci within that cluster. SPM j7~| denotes the f values for the 
stereotaxic coordinates. Areas shown in bold had significantly greater flow reductions with muscle afferent stimulation than with cutaneous 
stimulation. Stereotaxic coordinates of peak activation are expressed in millimeters and given in MNI space. Localization was based upon 
transformed Talairach coordinates (see text). BA, Brodmann's area. 
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Z = 44 



Z = 60 




Figure 3. Activations in response to cutaneous stimulation. Group brain activation in response to stimulation of cutaneous afferents of tine 
right index finger at 2 and 4T, using random effects analysis. Statistical images were coregistered with the SPM-MNl single subject T1 image. 
The axial/horizontal sections in the upper row show frontal and parietal activation areas. Sections in the middle row show areas of activation 
including the second sensory area (Sll). Sections in the lower row show an area of activation within both the right cerebellum. Colored bars 
labels indicate levels of activation in t value. For illustration purposes a cutoff of t = 3.5 has been used. 



Muscle afferent activation evoked significantly greater 
inhibition than cutaneous stimulation at multiple sites 
within the parietal lobes, parahippocampal gyri, temporal, 
occipital lobes, and one site in the cerebellum (Table 2). 
Cutaneous stimulation caused significantly more inhibi- 
tion at a single site in the contralateral parahippocampal 
gyrus (Table 4). 



4T versus 2T: muscle and cutaneous 

Muscle 4T stimulation evoked significantly greater activa- 
tion (f = 3.51-8.0) than 2T stimulation after masking, for 
five of the significant six sites in the ipsilateral cerebellum, 
the sites in the contralateral thalamus, globus pallidus, 
claustrum, insula, and one more superior site, in the 
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Afferent Activation of Sensorimotor Cortex and Cerebellum 



Table 3. All cutaneous stimulation excitation (random effects). Brain areas activated in response to supratlnreshold stimulation of the right 
index cutaneous afferents (f-test, P < 0.001 Uncorrected). 



Contralateral (Left) 



Ipsilateral (Riglnt) 



MNI coordinates 



Structure/Anatomy 



BA 



SPM .|7"| 



No. of 

voxels BA 



MNI coordinates 



No. of 

SPM {T\ voxels 



Frontal lobes 














Superior frontal gyrus 


6 


10 


12 


56 


4.5 


17 


Middle frontal gyrus 


8 


40 


20 


52 


6.23 


38* 


Precentral gyrus 


9 


34 


20 


44 


4.6 


38* 


Superior frontal gyrus 


8 


6 


30 


48 


4.96 


26* 


Medial frontal gyrus 


8 


6 


38 


44 


4.1 1 


26* 


Precentral gyrus 














Postcentral gyrus^ 


3 


62 


-22 


44 


4.07 


7 


Inferior frontal gyrus 














Superior frontal gyrus 


10 


20 


54 


16 


5.68 


26* 


Superior frontal gyrus 


10 


12 


62 


20 


4.24 


26* 


Inferior frontal gyrus 


46 


58 


38 


8 


4.66 


13* 


Inferior frontal gyrus 


47 


58 


36 


0 


3.9 


13* 


Parietal lobes 














Superior parietal lobule 


7 


12 


-78 


60 


4.58 


6 


Supramarginal gyrus 


40 


46 


-52 


32 


5.05 


6 


Temporal lobes 














Middle temporal gyrus 


39 


36 


-70 


28 


4.7 


18* 


Superior temporal gyrus 


39 


40 


-62 


24 


3.97 


18* 


Transverse temporal gyrus 


42 


62 


-16 


16 


5.19 


67* 


Middle temporal gyrus 














Middle temporal gyrus 














Middle temporal gyrus 














Insula 














Insula 


13 


48 


-18 


16 


5.47 


67* 


Insula 


13 


40 


6 


16 


4.6 


11 


Limbic lobe 














Parahippocampal gyrus 


28 


16 


-4 


-12 


4.83 


10 



Cerebellum 



9 38 8 44 4.5 18 

44 -46 10 20 4.14 17 



39 


-38 


-74 


16 


4.83 


56* 


39 


-58 


-68 


8 


4.58 


56* 


37 


-60 


-66 


0 


4.5 


56* 



-36 



-60 



-56 



Random effects model, second-level analysis. P< 0.001, minimum cluster size = 5. Areas shown in bold were significantly more activated by 
cutaneous stimuli than with muscle nen/e stimulation. Stereotactic coordinates of peak activations are in millimeters and given in MNI space. 
BA, Brodmann's area. 

'This point was initially allocated to area 4 by the daemon. 



postcentral cortex (Table 1). Sites within the ipsilateral 
parietal lobe, bilateral temporal lobes, contralateral fusi- 
form gyrus, and ipsilateral cerebellum were significantly 
more inhibited with 4T compared to 2T stimulation. 

No points were significantly more activated or inhib- 
ited by 4T cutaneous stimulation than with 2T within the 
significant regions of cutaneous activation. 

Discussion 

In this study, motor threshold was used as a substitute 
for muscle afferent threshold in contrast to the perceptual 
threshold used for cutaneous stimulation. Assuming that 



electrical threshold is inversely related to fiber diameter 
(Erlanger and Gasser 1937), it is likely that the 2T 
electrical stimulus activated most of the Group I muscle 
afferent fibers and A/? fibers. 4T stimulation would have 
activated most of the Group II and some Group III fibers, 
as well as the A8 fibers. The stimulus also activated motor 
axons and the ensuing muscle contraction would both 
directly and indirectly evoke muscle afferent volleys, sup- 
plementing the direct effects. Muscle afferents of Group I 
and II transduce muscle length and force and have 
important roles in proprioception (e.g., Proske 2006; 
Proske and Gandevia 2012) and the regulation of muscle 
activation both segmentally and supraspinally. Activation 
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Table 4. All cutaneous stimulation. Brain areas with reduced flow in response to stimulation of the right index cutaneous afferents (f-test, 
P < 0.001 Uncorrected, Voxels >5). 



Contralateral (Left) 



Ipsilateral (Right) 



MNI coordiates 



MNI coordinates 



Structure/Anatomy 



BA 



SPM |7| 



No. of voxels 



BA 



No. of 

SPM {Tl voxels 



Temporal lobe 

Superior temporal gyrus 
Limbic Lobe 

Parahippocampal gyrus 

Parahippocampal gyrus 
Subcortical 

Thalamus 



38 



35 



46 



10 



12 



32 



-32 



-20 -5.56 
16 -4.95 

-8 -4.12 



27 
83* 

83* 



30 



-38 



-4.12 



83* 



Combined data for both cutaneous nerve stimulation intensities showing reductions. Threshold for statistical significance was P< 0.001 
(uncorrected) with a cluster size of 5 or more voxels, using second level random effects analysis. Number of voxels indicates the overall spatial 
extent of significant clusters within the activated volume. Cluster sizes with asterisk (*) are foci within a larger cluster. SPM {7) denotes the 
f values for the stereotaxic coordinates. The area shown in bold had a significantly greater flow reduction with cutaneous stimulation than 
with muscle nerve stimulation. Stereotaxic coordinates of peak activation are expressed in millimeters and given in MNI space. Localization 
was based upon transformed Talairach coordinates (see text). BA, Brodmann's area. 



of cutaneous afferents, at 2 and 4 times perceptual thresh- 
old, was achieved by using digital nerve stimulation 
although some joint receptors were probably included. 
The largest cutaneous afferents mediate sensations of 
touch and vibration (Mountcastle 1974a). Projections of 
the large diameter primary afferents ascend to the thala- 
mus via the dorsal column-medial lemniscus and other 
pathways, the cerebellum via (for the upper limbs) the 
cuneocerebellar and rostral spinocerebellar and other 
tracts and brainstem targets (Mclntyre 1974; Sengul and 
Watson 2012). 

Selective electrical stimulation of both cutaneous and 
muscle afferent groups produced statistically significant 
increases in regional blood flow within cortical, 
subcortical, and cerebellar regions. Given we used single 
stimuli, likely to excite large diameter afferents, the areas 
of activation would be expected to coincide with those 
activated by cutaneous or muscle afferents with strong 
connectivity via oligosynaptic pathways. Responses 
corresponding to these projections might also be expected 
to be elicited by cutaneous stimulation or imposed move- 
ment under experimental conditions in intact animals or 
humans. The patterns of activation evoked by stimulating 
the two types of afferents differed. The cortical areas acti- 
vated with muscle stimulation included many associated 
with the generation of movement, including primary 
motor cortex (MI), thalamus, basal ganglia and cerebel- 
lum, as well as the insula, and other areas of the frontal 
and temporal lobes. Cutaneous stimulation activated mul- 
tiple areas in frontal, parietal, and temporal lobes, insula, 
parahippocampal gyrus, and cerebellum. 



Sensorimotor cortex 

Three significant foci of activation were present in the 
contralateral sensorimotor cortex following muscle affer- 
ent stimulation, with at least two distinct areas. MI lies in 
the precentral gyrus (area 4) and receives inputs from 
multiple areas including the premotor and supplementary 
motor area (SMA) (area 6). The activation evoked by pas- 
sive proprioceptive stimulation has been found to involve 
fewer cortical motor areas than active movement (Weiller 
et al. 1996; Mima et al. 1999) and in this study we, like 
these authors, found no areas of activation within nonpri- 
mary motor cortex such as SMA and premotor cortex. 
Brinkman and Porter (1979) reported that responses to 
passive movement were infrequent in SMA, at longer 
latency and insecure compared to those to the motor cor- 
tex. Naito et al. (2005) did find activation of nonprimary 
motor areas using muscle vibration, a process they sug- 
gested was due to muscle spindle afferent information 
being transmitted from "primary" motor and sensory 
areas to "non-primary" ones. This may have been facili- 
tated by the longer periods of continuous stimulation 
used (32s), allowing temporal summation of effects. 

Cortical activation following mixed nerve stimulation 
has been reported previously. Ibafiez et al. (1995) showed 
single focus in SI following median nerve stimulation at a 
similar height to our lower pair of foci. Although the pro- 
jection to area 3a is the densest cortical representation of 
muscle afferents, short latency projections to MI have 
been also been described (Colebatch et al. 1990). The MI 
representation for first dorsal interosseous muscle lies 
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within the hand area of the precentral gyrus and this can 
be locaHzed by its characteristic appearance of a "knob" 
(Yousry et al. 1997). Spiegel et al. (1999) used a sensory 
evoked potential paradigm to stimulate the median nerve 
at motor threshold and produced two distinct contralat- 
eral foci, lying in MI and SI. The center of activated cor- 
tex in MI lay 5 mm superior to the activated cortex in SI. 
Our most superior focus lay within the hand "knob" area 
for the right hemisphere of the image ch2 (Montreal 
Neurological Institute) and thus it is likely to represent 
activation within Ml by muscle afferents whereas the 
lower foci are likely to lie within SI. Our separation is 
greater than that reported by Spiegel et al. (1999) but we 
have confirmed that at least two separate sites of muscle 
afferent projections are present within the sensorimotor 
cortex, the superior one of which corresponds to the 
hand area of Ml as defined anatomically. 

We found that the Ml activation is significantly more 
powerful for muscle afferents than for cutaneous afferents. 
One area of significant activation associated with 
cutaneous stimulation lay within SI at the crown of the 
postcentral gyrus and thus in area 3b. Electrophysiological 
data suggest that the subdivisions of SI have different 
functional roles. Neurons in areas 3a and 2 respond to 
stimulation of deep receptors (Group I & II muscle 
afferents), whereas neurons in areas 3b and 1 exhibit a 
predominant responsiveness to stimulation of cutaneous 
receptors, consistent with our findings (Mountcastle 
1974b; Allison et al. 1991; Iwamura et al. 1993). 

Areas activated with muscle stimulation 
only 

Muscle afferent stimulation was associated with 
activations in the thalamus and basal ganglia. The two 
areas mapped as the claustrum may represent extensions 
of the bilateral insular activations and will not be consid- 
ered further. Muscle afferents from the upper limbs 
ascend in the dorsal columns, synapse in the cuneate 
nucleus and then terminate in the ventroposterior lateral 
(VPL) nucleus of the thalamus (Brodal 1981), consistent 
with the activation we found. The major projection from 
this nucleus is to postcentral sensory cortex - areas 1, 2, 
3 and 5 (Macchi and Jones 1997; Mai and Forutan 
2012). Parts of the "motor" thalamus, which projects to 
precentral motor areas, contain cells which respond to 
passive movements and muscle stretch at short latency 
(Vitek et al. 1994). Recordings from Vim in awake 
humans also demonstrate responses to active and passive 
joint movements, with neurons having characteristics 
consistent with muscle afferent input (Ohye et al. 1989). 

The contralateral globus pallidus showed activation 
with muscle afferent stimulation. The pallidum projects 
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to the VLa (VLo) area of the thalamus (Macchi and Jones 
1997) and labeling studies have shown projections from 
these specific thalamic nuclei to the striatum (McFarland 
and Haber 2001). In the somatosensory part of the GPi, 
up to 25% of neurones respond to kinesthetic stimuli in 
recordings made from patients operated on for 
Parkinson's disease (Lozano and Hutchison 2002). In 
monkeys, 22 and 37% of movement-related neurones in 
GPe and GPi respond to passive joint movement, but 
very few respond to cutaneous stimuli (DeLong et al. 
1985). 

Areas activated by cutaneous stimulation 
only 

Cutaneous stimuli, unlike muscle afferent activation, were 
associated with activation within the contralateral 
parahippocampal gyrus, bilateral temporal foci, contralat- 
eral parietal lobe, and multiple frontal areas, including 
ipsilateral sites. A number of studies have shown that 
cutaneous electrical stimulation can activate cortical areas 
in addition to the primary sensory cortex, including the 
contralateral superior parietal lobule (area 7), SMA/CMA, 
and insula as well as bilateral inferior parietal lobule 
(areas 39 and 40: Forss et al. 1994a; Ruben et al. 2001; 
Deuchert et al. 2002). The superior and inferior parietal 
lobules are considered to be multimodal association cor- 
tex and show widespread connectivity (Caspars et al. 
2012). 

Cerebellum 

The cerebellum was strongly activated ipsilaterally in this 
study following muscle afferent stimulation. High thresh- 
old muscle afferent activation was significantly more 
effective than low-intensity stimulation for nearly all the 
ipsilateral foci, supporting a role for higher threshold 
afferents. The cerebellum has a role in motor learning 
(Jenkins et al. 1994) and the ongoing modulation of 
motor output including timing of muscle contractions 
(Holmes 1952). Cerebellar activation has been shown 
previously using electrical stimulation of cutaneous and 
mixed nerves (e.g., Eccles et al. 1968). Painful 
intramuscular electrical stimulation of the brachioradialis 
muscle can produce ipsilateral activation in the declive of 
cerebellum, whereas painful cutaneous laser stimulation 
activates ipsilateral culmen (Svensson et al. 1997). Backes 
et al. (2000) showed significant activation in SI and the 
ipsilateral cerebellum with maximum current stimulation 
of the median nerve. Takanashi et al. (2003) reported cer- 
ebellar activation following stimulation of cutaneous 
afferents in the finger or toe. Muscle and other afferents 
project to the cerebellum through the dorsal and ventral 
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spinocerebellar tracts and their forelimb equivalents 
mainly to the intermediate and midline divisions of the 
cortex (Brodal 1981). Group la, lb, and Group II muscle 
afferent fibers project to neurons within both the dorsal 
spinocerebellar tract (DSCT) and ventral spinocerebellar 
tracts (VSCT) while the latter also receives input from 
high threshold (II and III) afferents of muscle, joint and 
cutaneous origin. While the DSCT terminates ipsilateraUy, 
afferents within the VSCT ascend in the contralateral cord 
and terminate bilaterally (Mclntyre 1974). Like Grodd 
et al. (2001) who studied voluntary movement, we found 
unilateral activation with both muscle and cutaneous 
stimuli and, like their report, there were at least two dis- 
tinct foci of activation within the ipsilateral cerebellum - 
one within the vermis and one more lateral. The latter 
was more difficult to localize and lay at the lateral bound- 
ary of the dentate nucleus, at the base of lobule VI 
(Dimitrova et al. 2002, 2006; Schmahmann et al. 1999). 

Short latency inputs from muscle and cutaneous affer- 
ents excite cerebellar neurons (Thach 1967; Eccles et al. 
1968, 1971). Muscle afferents caused significantly greater 
activation than cutaneous afferents in this study. This is 
consistent with the cerebellum's role in movement but 
contrasts with some previous reports using direct neuro- 
nal recordings from a similar region to our lateral activa- 
tion zone (paravermal lobules V and VI), where responses 
to probable cutaneous afferents (tapping) are more easily 
demonstrated using natural stimulation (Harvey et al. 
1977). These observations do not exclude muscle afferents 
arising from Golgi tendon organs (losif et al. 1972; Ishik- 
awa et al. 1972) or high threshold muscle afferents, both 
of which are unlikely to be activated by passive move- 
ments. Our findings indicate that peripheral afferent dis- 
charges are likely to contribute significantly to cerebellar 
activation occurring with voluntary movement. 

Insula 

The insula and surrounding somatosensory areas have 
been previously included as part of SII, the secondary 
somatosensory area, but this has now been subdivided 
(Kaas 2012). Activation within the insula region - the 
cortex lying in the floor of the lateral sulcus, occurred 
with both cutaneous and muscle afferent stimulation. It 
was contralateral only for cutaneous stimulation but bilat- 
eral for muscle afferent stimulation. It was significantly 
greater following muscle than cutaneous afferent stimula- 
tion, and significantly greater for 4T than 2T stimulation 
implying a role for Group II and probably Group III 
afferents. The posterior insula is a multisensory area 
(Mazzola et al. 2006) and stimulation can evoke painful 
sensations (Ostrowsky et al. 2002; Afif et al. 2008). It may 
code the magnitude of sensory inputs (Baliki et al. 2009). 
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We did not match our stimuli for their subjective dis- 
comfort and the 4T muscle stimulation was considered to 
be the most painful. 

Previous imaging studies have looked at the involve- 
ment of muscle and cutaneous afferents in pain. Svensson 
et al. (1997) used high-frequency intramuscular electrical 
stimulation (brachioradialis m.) and cutaneous stimula- 
tion with PET imaging to observe regional activation in a 
variety of mainly contralateral cortical regions which 
overlapped for the two modalities. They reported that 
rCBF in the anterior insular cortex increased with nox- 
ious cutaneous stimulation. Electrical stimulation of cuta- 
neous digital afferents may activate contralateral SI and 
SII areas with somatotopic representation (Deuchert et al. 
2002). A study by Niddam et al. (2002) compared painful 
to non-painful intramuscular electrical stimulation using 
event-related fMRI (3T), and found that painful muscle 
stimuli increased activation in ipsilateral SII and IPL 
(areas 43 & 40), middle and posterior insula as well as 
contralateral middle frontal gyrus (area 10), precentral 
gyrus (area 44), superior temporal gyrus (area 22), limbic 
lobe (areas 23, 24 and 32), anterior and posterior insula, 
and caudate body and tail. Activation within the posterior 
insular cortex was significantly greater for painful muscle 
stimulation. 

Temporal areas 

Areas within the middle and upper temporal lobes were 
activated during both muscle and cutaneous afferent 
stimulation. Area 42 and area 41 correspond to Heschl's 
gyrus, the primary auditory field (Brodal 1981). Areas 22, 
37 and 39, all showing activation, are perisylvian language 
areas. None of these areas is directly related to somato- 
sensory input, but they have shown activation in previous 
studies. Niddam et al. (2002) reported bilateral activation 
within area 22 (right predominant) in response to both 
non painful and painful muscle stimulation, an effect that 
attributed to stimulus "salience". In this study, for muscle 
afferent stimulation, there was a clear left-sided predomi- 
nance. As fMRI imaging is noisy and semirhythmical, it is 
possible that the increased activation in these areas repre- 
sented heightened awareness of external stimuli during 
the period of stimulation. 

Relationship with evoked potentials 

Ironically, given the wide range of sites of activation, early 
studies using evoked potentials showed responses in cere- 
bral cortex only when stimuli engaged Group II or Group 
III volleys, leading to the conclusion that Group I activity 
alone was ineffective (Mclntyre 1974). In contrast, affer- 
ent volleys from skin and joint nerves readily elicited 
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cortical potentials. Gandevia and Burke (1988, 1990), 
using motor point and intrafascicular stimulation, 
reported evoked potentials for both proximal and distal 
muscles of the arm, but the potentials were small com- 
pared to those for a digital nerve. Allison et al. (1991) 
concluded that short latency SEPs evoked by stimulation 
of the median, a mixed nerve, were primarily the result of 
cutaneous fibers and that muscle afferents contributed lit- 
tle. Our findings suggest that evoked potentials are not a 
reliable means by which to assess the overall cortical pro- 
jections of afferent species. The larger evoked responses to 
cutaneous afferents are likely to be due in part to the 
location and orientation of the generators, the relative 
number of cutaneous afferents and perhaps the synchrony 
of the voUey. Evoked potentials will be biased toward cor- 
tical responses whereas here we have compared muscle 
and cutaneous afferent effects at the cortical, subcortical, 
and cerebellar levels. 

The relationship between changes in BOLD signal, as 
recorded here, and electrophysiological measures has been 
the subject of a number of investigations. EEG potentials 
are believed to arise from summed postsynaptic potentials 
and horizontally oriented cortical cells, and subcortical 
structures contribute little, if anything, to the scalp- 
recorded EEG (Aminoff 1999). Likewise, simultaneous 
recordings in the visual cortex of anesthetized monkeys 
have shown that BOLD signal increases reflect increases 
in neural activity, which is more closely related to local 
field potentials (LFPs), a measure of afferent input, than 
to neuronal discharges (Logothetis et al. 2001). This may 
apply particularly to high-frequency LFPs (100-150 Hz, 
Ojemann et al. 2010) but may vary for different cortical 
regions (Huettel et al. 2004). Our findings indicate a 
widespread influence of muscle afferents, including in the 
sensory cortex, which is more intense than that for 
cutaneous afferents. The effects are not simply increases 
in areas of BOLD activation, as muscle afferents also 
cause significantly larger areas of reduced activation than 
cutaneous afferents. Cortico-cortical inhibition of nonsa- 
lient sensory receptive areas is a general feature of sensory 
processing and appears to be mediated by GABAergic 
transmission (lurilli et al. 2012). 

Functional significance 

Studies using active movement must go to considerable 
lengths to avoid differences in the kinematics of the 
movements so as not to confound their results by differ- 
ences in the peripheral afferent activations evoked by the 
movements themselves (e.g., Schaechter and Perdue 
2008). Disorders of sensorimotor integration are common 
in movement disorders (Tempel and Perlmutter 1990; 
Abbruzzese and Berardelli 2003). In dystonic disorders. 
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sensory input may either trigger or improve the 
abnormality and the widespread projections we have 
demonstrated, including the basal ganglia, could underlie 
these phenomena. Clinically, it is widely accepted that 
normal movement is not possible in the absence of pro- 
prioceptive sensation and deafferented patients typically 
perform ataxic movements ("sensory ataxia" - e.g.. 
Holmes 1952) and this may be explained by the wide- 
spread projections of muscle afferents to motor-related 
areas interacting with the motor command. Our findings 
also support Holmes' observation that it is not simply a 
loss of conscious proprioception that underlies the clinical 
deficits but also the concomitant loss of the subconscious 
afferent input to the cerebellum and other areas. Func- 
tional imaging indicates that the central effects of muscle 
afferents are more widespread and more focused on the 
motor system than are those of cutaneous afferents. The 
same applies to projections to sensorimotor cortex but is 
contrary to the evidence that might be deduced from 
evoked potentials. 
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